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Abstract 

Rotating transonic flows are long known to admit standing or oscillating shocks. The 
excess thermal energy in the post shock flow drives a part of the in falling matter 
as bipolar outflows. We compute mass loss from a viscous advective disc. We show 
that the mass outflow rate decreases with increasing viscosity of the accretion disc, 
since viscosity weakens the centrifugal barrier that generates the shock. We also 
show that the optical depth of the post-shock matter decreases due to mass loss 
which may soften the spectrum from such a mass losing disc. 
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1 Introduction 



Accretion onto black holes has been intensely studied for the last three decades, 
in order to explain observed luminosities of galactic black hole candidates and 
AGNs, their spectral states, and the formation of jets or outflows around black 
hole candidates. Since black holes have neither hard surface, nor intrinsic at- 
mosphere from which outflows may originate, jets /out flows has to originate 



from the accreting matter. Recently. iGallo et al. I ( 120031 ) showed that jets orig- 



inate from accretion discs which are in hard spectral states. This fact suggests 
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that the origin of jet is strongly correlated with the spectral state of the ac- 
cretion disc, and therefore with the accretion process itself. 



The unique inner boundary condition for accretion onto a black hole, neces- 
sarily demands that the black hole accretion is globally transonic. IChakrabarti 
( 119891 ) showed that transonic matter suffers centrifugal pressure mediated 
shock, and the shock is located further out as the specific energy and the 
specific angular momentum of the flow is increased. Such a shocked accretion 
disc m odel was used to compute the spectral states of the bla c k hol e candi- 
dates ( Chakrabarti &: Titarchuk . 19951 : Chakrabarti fc Mandal . 2006 ). where 
the post-shock tori — CENBOL (CENtrifugal Pressure supported Bound- 
ary Layer), produces the hard power-law tail by inverse- Comptonizing the 
softer photo ns from the outer disc. Mo reover, it has also been shown , both 



numericallv (IMolteni et al. I. Il994j. Il996l) an d analytically (jChakrabartil . Il999l : 



Das fc Chakrabarti! . 1 19991 : iDas et al. 1 . 1200 ll ) , that the extra thermal gradient 



force in the CENBOL drives a sign i ficant portion of accreting matter as bipolar 
outflows. In particular. iDas et al. I ( l200ll ) showed that outflows from CENBOL 
are generated in hard spectral states of accretion discs. These bipolar outflows 
are believed to be the precursor of astrophysical jets. 



So far, there is no theoretical attempt to compute mass outflow rates from 
a viscous advective disc. As viscous flow approaches the black hole, the local 
angular momentum of the flow decreases while the local energy increases. 
Thus, in a viscous flow there are two competing processes to determine the 
shock formation. Should the shock move inwards as we increase the viscosity 
for flows with identical outer boundary conditions or vice-versa? How would 
the position of shock affect the mass outflow rates? We are going to address 
the above mentioned issues in the present paper. 



In the next section, we present the model assumptions and the equations of 
motion of the disc-jet system. In §3, we present the solution procedure. In §4, 
we present the results, and in the last section we draw concluding remarks. 



2 Model Assumptions and Equations of motion 



We begin with a steady, viscous, axisymmetric accretion flow on to a Schwarzschild 
black hole. The space-time property is approximated by the Paczyhski-Wiita 



potential (jPaczyhski fc Wiita I . Il980l ). In addition, jets are considered to be 



inviscid. Outflows are supposed to be lighter than the accretion disc, and def- 
initely of lower angular momentum than at least the outer edge of the disc. 
Consequently, the differential rotation is expected to be much less than that 
in the accretion disc. Therefore, viscosity in jets can be ignored. In Fig. 1, 
a schematic diagram of an advective accretion disc-jet system is presented. 
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Fig. 1. Schematic diagram of disc-jet system (with up-down symmetry) where the 
position of black hole (BH), shock location {xg), CENBOL, centrifugal barrier (CB), 
funnel wall (FW) etc are shown. The dotted line is the axis of symmetry. 

In this figure the pre-shock flow, the shock location (xs), the post-shock flow 
(abbreviated as CENBOL) and the jet between the funnel wall (FW) and 
centrifugal barrier (CB) have been marked. We define FW and CB in §2.2. 
Furthermore, we use a system of units such that 2G = Mbh = c = 1, where 
G, Mbh, and c are the gravitational constant, the mass of the black hole, and 
the velocity of hght, respectively. 



2. 1 Equations for accretion 



In the steady state, the dimensio nless hydrodynamic equations of motion for 
accretion are (jChakrabartil . Il996l ). 



the radial momentum equation : 



du IdP X^(x) 
u— + + 



dx ' p dx x^ 2{x — 1)^ 
the baryon number conservation equation : 

M = 27rSMx, 
the angular momentum conservation equation : 



u- 



dXjx) 
dx 



(la) 



(16) 



(Ic) 
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Fig. 2. Jet geometry for A = 1.75. In figure OA = rpw^ OC = rcB, QC = xcb, 
CM = ycB] CM = BP{yj) = AN{yFw)] QA = xpw; QB = xj; And OB = rj. 
FW and CB are marked in the figure. 

and finally, the entropy generation equation : 

uT- = --x{ga + m {Id) 

where, variables m, a, p, P and and \{x) in the above equations are 

the radial velocity, sound speed, density, isotropic pressure, angular velocity 
and specific angular momentum of the flow respectively. Here S is vertically 
integrated density, M4(/)(= — anH, 11 is the vertically integrated total pressure, 
an is the viscosity parameter) denote the viscous stress, s is the speciflc en- 
tropy of the flow, and T is the local temperature. The disk is assumed to be 
in hydrostatic equilibrium in the vertical di rection. Furthermore, q = In+i/In, 



n = 1/(7 - 1) and 4 = (2"n!)V(2n + 1)! flMatsumoto et al. I . Il984f ). and the 



adiabatic index 7 = 4/3 is used through out the paper. 



2.2 Equations of motion for outflows 



Numerical simulations by iMolteni et al. I (119961 ) suggests that the outflowing 



matter tends to emerge out between two surfaces namely, the funnel wall (FW) 
and the centrifugal barrier (CB). In Fig. 2, the schematic diagram of the jet 
geometry is shown. Geometric surfaces CB and FW are marked in the flgure. 
The centrifugal barrier (CB) surface is deflned as the pressure maxima surface 
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and is expressed as 

xcB = [2A^rcB(rcB - 1)J , (2a) 

where, tcb = \/xcb + ycB= spherical radius of CB. Here xcb, Vcb are the 
cyhndrical radius and axial coordinate {i.e., height at tcb) of CB. We compute 

the jet geometry with respect to ucb i-^., Ufw = Vj = Vcb, where ypw and 
■ijj are the height of FW and the jet at tcb, respectively. The FW is obtained 
by considering null effective potential and is given by 



. (A^-2) + V(P-2)^-4(l-y^,) 

where, xpw is the cylindrical radius of FW. Using Eqs. (2a) and (2b), we 
estimate the cylindrical radius of the outflow 

Xpw + XcB /P, X 

= ^ ■ (2c) 



The spherical radius of the jet is given by rj = yxj + yj. In Fig. 2, OB (= Tj) 
defines the streamline (solid) of the outflow. The total area function of the jet 
is obtained as, 

A = 27r(a;^B - x^y^) {2d) 
The integrated momentum balance equation for jet is given by, 

1 1 
= + + ______ (3a) 

where Sj and Xj are the specific energy and the angular momentum of the jet, 
respectively. The integrated continuity equation is, 

Mont = PjVjA, (36) 

and instead of the entropy generation equation we have the polytropic equation 
{pj = KjpJ) of state for the jet. In Eqs. (3a-3b), the suffix 'j' indicates jet 
variables, where Vj, aj, and pj are the velocity, the sound speed and the density 
of the jet. 



3 Method 



Since post-shock disc is the source of outflows, we are interested in only those 
accretion solutions which includes stationary shocks. A shocked accretion so- 
lution has two X-type critical points — inner (xd) and outer (xco) critical 
points. Critical points are radial distances at which the gradient of veloc- 
ity i.e., du/dx^O/0. The accretion solution is solved by supplying Xd and 
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angular momentum Aj at Xd and integratinK the equations once from Xd in- 
wards and then outwards [as in (jChakrabarti fc Dad . |200J)]. The conditions 
for steady shock require conserva tion of energy flux, rnass flu x, and momentum 
flux across the shock front [e.g., ( iLandau fc Liftshitzl . Il959l )] and are given by, 



M 



0; [n] = 0. 



(4a) 



In presence of mass loss, the condition for conservation of mass flux takes the 
following form, 

M+ = M_ - Mo,t = M-(l - Rrn), m 

where, subscripts "— " and "+" refer to quantities before and after the shock, 
respectively. The mass outflow rate R^n is the ratio between mass flux of the 
outflow (Mout) and the pre-shock accretion rate (M_). Using the above shock 
conditions, the pre-shock sound speed and bulk velocity can be expressed as. 



7m: 



(5a) 



and. 



1 - 



27 



2Ci 



Xq 



(56) 



where Ci = (1 — Rm){ga'^ + 7f^)/M+, and {Ss, K) are the energy and the an- 
gular momentum at the shock. As we integrate from Xd outwards, Eqs. (5a-b) 
are used to compute pre-shock flow variables. These pre-shock variables [i.e., 
U-, a_) are then used to find Xco- However, Eqs. (5a-b) show that the infor- 
mation of mass loss is in the shock condition itself. This implies that, in order 
to obtain shock the accretion-jet equations have to be solved simultaneously. 



Jet equations are s olved by employing the so-called 'critical point analysis' 
( jChakrabarti 119901 ). We differentiate Eqs. (3a-3b) with respect to r(= vcb), 
and obtain the expression for dvj/dr. The critical point (r^c) condition for jet 
is. 



1 



'drj_ 
dr 



/ dxj 
dr 



X 




-1 



(3c) 



Once the flow variables at the jet critical points are known, we integrate the 
jet equations from Vc to the jet base (xs). 

The expression for R^ is calculated by assuming that jets are launched with 
the same density as the post-shock flow and is given by. 



Rn 



Rvj{xs)A{xs 



M_ 



AttJ-x 

V 7 



2" 3/2, 



— l)a+'U_ 



(6) 
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Fig. 3. Upper panel: Inflow Mach number (M = u/a) with log(x). The inflow param- 
eters are Xd = 2.445, \i = 1.75, and an = 0.002, where, Xg = 19.21, Sg = 0.0014, 
Xs = 1.756 Xco — 206.49, Xq — 1.772. The dotted curve is the shock free solution. 
Lower panel: Outflow Mach number (Mj = Vj/aj) with log{xj), the outflow critical 
point Xjc = 81.26 {rjc = 346.18), and the jet coordinates at the base is given by 
Xjb = 10.27 {rjb = 16.62). The mass loss rate is Rm = 0.055. 

where the compression ratio is expressed as R = E_|_/E_. The self-consistent 
method to obtain shocks in presence of mass loss is the following : initially 
the mass loss is not considered, and the shock is found out for accretion flows. 
Using {£j, Xj) = {Eg, Xg) as inputs, we solve the jet equations and calculate the 
values of jet variables at r = Xg. Once we obtain the value of Vj, A at r = Xg, 
we use Eq. (6) to calculate i?^. The value of Rm is used in Eqs. (5a-5b), and 
the new shock location is obtained. We iterate this process till the solutions 
converge. 



4 Results 

Supersonic matter, which was subsonic at the outer edge of the disc, suffers a 
discontinuous transition due to centrifugal barrier at shock {xg) and becomes 
sub-sonic. A significant part of the post-shock matter is ejected as bipolar 
outflows due to excess thermal pressure at the shock and the rest enters into 
the black hole through Xd- In Fig. 3, we present such a global inflow-outflow 
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Fig. 4. Variation with Si^, for an = (solid), 0.005 (dotted), 0.01 (dashed), and 
0.015 (dashed-dotted) , 0.02 (fong-dashed) respectively, and Aj = 1.75. 

solution. In the top panel the Mach number of the accretion flow is plotted 
with log{x). The solid curve represents shock induced accretion solution. The 
inflow parameters are Xd — 2.445 , Aj = 1.75, and an = 0.002. In the lower 
panel, the outflow Mach number Mj is plotted with log{xj). In presence of 
mass loss {Rm = 0.055) the shock forms at Xg = 19.21 (vertical line in the top 
panel), and the outflow is launched with energy and angular momentum at 
the shock (£s, K = 0.0014, 1.756). The outflow is plotted up to its sonic point 
{xjc^ S1.26). 

We now present how the mass outflow rate is affected by viscosity in the disc. In 
Fig. 4, is plotted with (the energy at Xd), for a set of ctn- For instance, 
solid curve (an = 0) represents mass outflow rates for inviscid accretion. 
Other curves are obtained for the following values of an: 0.005 (dotted), 0.01 
(dashed), and 0.015 (dashed-dotted), and 0.02 (long dash-dotted) respectively. 
All the curves are drawn for Aj = 1.75. We observe that higher energetic flow 
produces higher i?^ for a given an- Conversely, for the same value of if 
an is increased, then i?^ is decreased signiflcantly. It is to be noted that, for 
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each an there is a cut-off in R,n at the higher energy range, since standing 
shock conditions are not satisfied there. Non-steady shocks may still form in 
those regions, but investigation of such phenomena is beyond the scope of this 
paper. 

In Fig. 5a, i?™ is plotted with an, for {xd, Xi) = (2.313, 1.8) (solid), (2.375, 1.775) 
(dashed) and (2.445, 1.75) (dotted), respectively. For a given an, Rm increases 
with higher Aj. This suggests that outflows are primarily centrifugal pressure 
driven. However, i?^ decreases drastically with an for each of the curves. The 
question however is, why decreases with an? We wish to discuss the under- 
lying physics and consider a representative case [e.g., {xd, Xi) = (2.445, 1.75)]. 
In Fig. 5b, shock location (xg) is plotted with an, for solutions without mass 
loss (dotted) and with mass loss (solid). As viscosity is increased the shock 
location is decreased. Generally, in presence of mass loss shock forms closer to 
the black hole. Mass loss reduces the post-shock pressure, and therefore the 
shock moves inwards to maintain pressure balance across it. However, the two 
curves tend to merge at higher an, which signifies that the mass loss decreases 
with an- In Figs. 5(c-d), and Sg are plotted with an, respectively. In this 
case, Xi is same, therefore with increasing an there is a moderate increase of 
X{x) as we integrate outwards, however, S{x) decreases sharply, which results 
in the lower value of Xg. The fractional increase of A^ is ~1.1%. The frac- 
tional decrease of Ss is ~8.2%. As Xg moves inwards, the base area of 
the jet decreases, and consequently the base velocity [vj\r=xs) of the jet also 
decreases, while the compression ratio (R) increases. In Fig. 5e, VjA (solid) 
and it! (dotted) is plotted with an- The quantity VjA decreases by 90%, while 
R increases by mere 60%. From Eq. (6), we see that in the denominator of 
the expression for Rm, there is a term xl^'^{xg — l)a_|_-u_. In Fig. 5f, both the 
quantities x^^'^^Xg — l) (dotted) and a^u^ (solid) is plotted with an- The quan- 
tity x'l^'^^Xg — 1) decreases by 37% and a+M_ (multiplied by a factor 5x10^ 
in the figure) increases by 42%. Therefore, as Xg decreases, the decrease of 
VjA dominates all other quantities in Eq. (6), and hence Rm decreases with 
the increasing an- In Figs. 5(b-f), we have fixed {Xi,Xci), and increased an to 
see how viscosity affects shock and consequently Rm. If one increases an for 
fixed {Xi,Xci), then will decrease anyway [e.g., Fig. (Ic) of Chakrabarti & 
Das 2004], hence Eg will be lower for higher an- In other words, we have stud- 
ied less energetic accretion flows though we have increased the an parameter, 
and saw that for such flows increasing viscosity parameter decreases the shock 
location and therefore the mass outflow rate. The question however is, what 
happens to flows starting with same outer boundary conditions, as viscosity 
is increased? 

In Fig. 6, the variation of M (upper panel), A (middle panel), and S (lower 
panel) with log{x) is presented. The flow is launched at the outer edge (xinj = 
200) with energy and angular momentum corresponding to (^inj, Ainj) = (0.002, 1. 
For inviscid case, the shock forms at Xg — 20.3 (dashed) and mass outflow 
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Fig. 5. (a) Rm vs an, for for Aj = 1.8 (solid), 1.775 (dashed) and 1.75 (dotted), 
respectively. Inner sonic points are Xd = 2.313, 2.375, and 2.445, respectively, (b) 
Variation of Xg with ccn, \ = 1-75 Xd = 2.445. Dotted curve represents solution 
without mass loss, and solid represents Xg with mass loss, (c) Variation Xg with 
«n,Ai = 1.75 Xci = 2.445. (d) Variation 8s with au,K = 1-75 Xd = 2.445. (e) 
Variation of R (dotted) and vjA (solid) with an, Aj = 1.75, Xd = 2.445. (f) Variation 

of x^^^{xs — 1) (dotted) and as+u-x5e + 4 (solid) with an, Aj = 1.75, Xd = 2.445. 



10 



1.7 
0.003 



EiD 0.0025 



\ 

" \ 






" i : I"- 
- i i 1 

X i 1 


"\. 




\. i 1 

\ V ! 


^\ 




: \. \ : 




















I.I. 


1 





log(x) 



Fig. 6. Variation of M (upper panel), A (middle panel), £ (lower panel) with log{x), 
for an = (dashed) and an = 0.003 (dotted). The outer edge quantities are 
Xinj = 200, (£^inj,Ainj) = (0.002,1.75). The shocks are Xg = 20.3 (dashed) and 
Xs = 11.65 (dotted). 

rate is calculated to be Rm = 0.093. As viscosity is incorporated (an = 0.003; 
dotted), the shock forms at Xs = 11.65 and the outflow rate is reduced to 
Rm = 0.07. As viscosity is enhanced, S{x) is increased and \{x) is decreased 
along the flow. At the shock, the effect of energy increment, cannot compensate 
the effect due to angular momentum reduction. This weakens the centrifugal 
barrier and causes the shock to move inward. As shock moves inwards the 
jet-base area decreases, which reduces the mass outflow rate. 

So far, we have shown that the mass outflow rates strongly depends on shock 
location Xg, and formation of shock itself depends crucially on the viscosity 
parameter. We have also shown that shocks form closer to the black hole for 
viscous flows than in inviscid accretion flow. However, we are yet to comment 
on the observational consequence of such phenomena. Admittedly, it is beyond 
the scope of the present paper to compute the detailed spectrum of the flow, 
but we would like to make some qualitative remarks on the spectral proper- 



ties of such flows. It is to be remembered, IChakrabarti fc TitarchukI (1l995l ) 
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Fig. 7. Variation of r with (x), for {an,Rrh) = (0,0) (solid), (an,^m) = (0,0.093) 
(dotted), ian,Rm) = (0.003,0) (dashed), and (an,-Rm) = (0.003,0.0.07) (dashed- 
dotted). The outer edge quantities are xinj = 200, (finj, Ainj) = (0.002, 1.75). Here r 
is plotted for CENBOL. 

showed that hard power law photons were generated at CENBOL by inverse- 
Comptonization of soft photons from the pre-shock region of the disc. Thus 
we calculate the optical depth of the CENBOL, as this provides a qualitative 
understanding of the nature of the spectrum. The definition of the optical 
depth for photons entering the CENBOL is given by, 

/X 
Kpdx', 
Xs 

where, k is the opacity assuming Thompson scattering cross-section, and p 
is the density of the flow. We consider matter with identical outer boundary 
conditions (xinj, £"111], Ainj) = (200,0.002,1.75). The mass accretion rate at Xinj 
is M_ = 0.3MEdd; and the central mass is Mbh = IOMq. In Fig. 7, t{x) is 
plotted only for the post-shock disc with {au,Rm) = (0,0) (solid), (0,0.093) 
(dotted), (0.003,0) (dashed), and (0.003,0.07) (dashed-dotted), respectively 
For the most simplified case (solid), the shock is at Xg = 25, and the optical 
depth T~l is at Svg. Therefore most of the soft photons entering the CENBOL 
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fails to sample the inner part of the disc. In presence of mass loss (dotted), 
the density of the CENBOL is lowered which reduces its optical depth. For 
viscous flow (dashed), optical depth is much lower compared to the inviscid 
flow (solid) as the shock forms {xs = 14.2) closer to the black hole. If we allow 
mass loss (dashed-dotted), the optical depth is reduced further. Therefore, 
most of the soft photons coming from the pre-shock region could interact with 
the hot electrons of the inner part of the disc (r ~ 1 at x ^ 4rg) and there by 
cool it. A cooler CENBOL will evidently produce softer spectrum compared 
to the inviscid case. Chakrabarti (1998) calculated spectrum of a inviscid disc 
which is loosing mass from the post-shock region, and indeed he found that the 
spectrum softens, which correspond to solutions similar to the dotted curve. 
From our self-consistent study, we predict that the twin effects of viscosity 
and mass loss will soften the spectrum even further. 



5 Concluding Remarks 



Early studies of transonic accretion flow showed the existence of standing or os- 
cillating shocks (IChakrabartl Il989l ) . The presence of shoc ks were also thought 

to be the real cause behind hard state of the disc spectrum (IChakrabarti fc Titarchukl . 
19951 ) ■ and also generating outflows (IChakrabartil . Il999l ). Studies of shocks in 
viscous accretion disc were also undertaken, but were confined to obtain the 
global solution for accretion fl ows with or without shocks and/or the detaile d 
study of the parameter space (IChakrabartil . Il996l : IChakrabarti fc Dad . |2004| ). 
Important as they are, since a proper understanding enhances the grasp on the 
physics of such flows, no attempts have been undertaken to link these viscous 
flows with observables such as outflows and sp ectrum. In this p aper we have 



extended our earlier computation of mass l oss (IDas et al. I. l200ll). bv emplov- 



ing the knowle dge of viscous dis c solutions (IChakrabarti fc Dad . 120041 ) . In the 



earlier paper (IDas et al. I . l200ll ). the disc was considered to be inviscid and 



the jet was non-rotating. Furthermore, the shock condition was not modified 
for mass loss. Presently, we have included viscosity and considered rotating 
outflows, in particular we have concentrated on studying how viscosity affects 
shocks and in turn, how shocks affect the mass outflow rates. 



We have shown that mass loss from the post-shock disc reduces post-shock 
total pressure, and hence the shock front moves towards the black hole. We 
have also shown that the mass outflow rate depends on the energy, as well 
as the angular momentum of the flow, and increasing both, in turn, increases 
the mass outflow rate. However, less energetic accretion will produce lower 
mass outflow rates, even though angular momentum is increased outward, by 
increasing the viscosity parameter. 



More interestingly, we have shown that outflow rate decreases with the increase 
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of viscosity parameter for flows with same outer boundary conditions. This 
is due to the weakening of centrifugal barrier with the increase of viscosity 
parameter. We have also shown that the post shock optical depth decreases in 
presence of viscosity and mass loss, which may soften the emitted spectrum. 
A detailed study of spectrum from such viscous discs, by including cooling 
processes is under consideration and will be reported elsewhere. 
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